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ABSTRACT 

The analysis of the light curves of 48 B-type stars observed by Kepler is presented. Among 
these are 15 pulsating stars, all of which show low frequencies characteristic of SPB stars. 
Seven of these stars also show a few weak, isolated high frequencies and they could be con- 
sidered as SPB//? Cep hybrids. In all cases the frequency spectra are quite different from what 
is seen from ground-based observations. We suggest that this is because most of the low fre- 
quencies are modes of high degree which are predicted to be unstable in models of mid-B 
stars. We find that there are non-pulsating stars within the /? Cep and SPB instability strips. 
Apart from the pulsating stars, we can identify stars with frequency groupings similar to what 
is seen in Be stars but which are not Be stars. The origin of the groupings is not clear, but may 
be related to rotation. We find periodic variations in other stars which we attribute to proximity 
effects in binary systems or possibly rotational modulation. We find no evidence for pulsating 
stars between the cool edge of the SPB and the hot edge of the S Set instability strips. None 
of the stars show the broad features which can be attributed to stochastically-excited modes 
as recently proposed. Among our sample of B stars are two chemically peculiar stars, one of 
which is a HgMn star showing rotational modulation in the light curve. 
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1 INTRODUCTION 

The Kepler mission is designed to detec t Earth-like planet s around 
solar-type stars by the transit method dKoch et al]|2010h . Kepler 
will continuously monitor the brightness of over 150000 stars for 
at least 3.5 yr in a 105 square degree fixed field of view. In the 
course of the mission, extremely accurate photometry will be ob- 
tained. The mean top-of-the-noise level in the periodogram (which 
we will call the grass level) for a star of magnitude 10 is about 
3 ppm for an observing run of 30 d, dropping to about 1 ppm after 
one year. For a stai^ of magnitude 12 the corresponding grass noise 
level is about 6 ppm for one month and 3 ppm after one year. These 
values apply to a frequency range 0-20 d~^, the noise level drop- 
ping slowly towards higher frequencies. With this unprecedented 
level of precision one might expect to discover new behaviour for 
many stars. 

Stability analyse s lICox et alj 

iDziembowski & Pamvatnvkhl 1 19931 : iDziembowski et al] 



1993 ; 



1993h 



show that an opacity bump due to a huge number of absorption 
metal lines can destabilize low-order p and g modes in /? Cephei 
stars and high-order g modes in slowly-pulsating B (SPB) stars. 
Driving of pulsations due to the k mechanism can only occur if 
certain criteria are met. One of these criteria is that the pulsation 
period is of the same order or shorter than the thermal timescale. 
Otherwise the driving region remains in thermal equilibrium and 
cannot absorb/release the heat required for driving the pulsations. 
Another requirement is that the pressure variation is large and 
varies only slowly within the driving region. This criterion is 
satisfied for low radial order p modes and (in some models) 
high-order g modes with I < 6 in /3 Cep stars in the iron-opacity 
bump region. In /3 Cep stars this opacity bump is located in 
a relatively shallow layer where the thermal timescale is well 
below 1 d. It turns out that the timescale constraint is fulfilled 
on ly for p modes of low radial order and for g modes with I > 
6 ( IDziembowski & PamvatnYklilll993l) . In less massive stars, the 
iron-opacity bump is located in a deeper layer and, in addition, the 
luminosity is lower. Both factors contribute to an increase in the 
thermal timescale. Consequently, all p modes are stable because 
their periods are much shorter than the thermal timescale, while 
high radial order g modes with low / fulfill the timescale and 
pressure amplitude constraints and can now be driven. Thus we 
find two instability regions among the B stars: the high-luminosity 
/3 Cep instability strip and the low-luminosity SPB instability strip. 

A number of /3 Cep stars, however, sh ow low-fr equency pul- 
sations characteristic of SPB stars (see e.g . , 
Chapellier et al.l 120061 lDe"Cat et al.l I2OCT : 



Le Contel et al. 200 li: 



Pigulski & Poimaiiskil 

2008al lbl: lHandleJl2009h . These are known as /3 Cep/SPB hybrids. 



In general, high radial order g modes are heavily damped in B stars 
of high mass because they have high amplitudes in the deep lay- 
ers just above the convective core. The temperature variation gra- 
dient in this region is very large, leading to significant heat loss 
and damping of the pulsations. However, the structure of the eigen- 
functions plays an important role and for some modes driving in 
the iron-opacity bump region exceeds damping in the inner layers. 
Hence some low-frequency g modes can occur together with high 
frequency p modes in the same star. 

Unfortunately, the g modes observed in some /3 Cep 
stars are usually found to be stable in the models (e.g. 
IDziembowski & Pamvatnvkl]|2008l) . It is also found that the ob- 
served range of instability at high frequencies is wider than can be 
explained. The problem may be resolved by postulating an over- 
abundance of iron-peak elements in the driving region by up to 



a factor of four jPamvatnvkh et al.ll2004l). The use of OP op aci- 
ties and the recently-revised solar mixture dAsplund et al.l2009l) re- 
quires a modest 50 percen t enhancement of Fe in the iron-opacity 
bump for u Eri models ( iZdravkov & Pamvatnyk h 2008). Since 
there is no observational evidence to suggest that the current metal 
abundances are in error by that amount, the required opacity en- 
hancement needs the increase in metal abundance to be confined 
to the iron-opacity bump region. Unfortunately, any enhancement 
of Fe in this regi on cannot be hidden and should be detectable in 
the photosphere ( ISeatonlll999l) . Either a further upwards revision 
of the opacities is required or some physics is missing from model 
calculations. 

The hybrid stars provide an important test of the models be- 
cause they probe a very wide frequency range. Hybrid behaviour is 
more common in m odels using OP opacities than OPAL opacities 
( iMigUo et a"l]|2007h . This is mainly due to the fact that the metal 
bump opacity in the OP data occurs at a deeper level where the 
temperature is about 20 000 K higher than in OPAL. Also, the pe- 
riod gap betwe en /3 Cep stars and SPB stars is expected to decrease 
as I increases jBalona & Dziembowskj||T999l) . If most of the low- 
frequency modes in /3 Cep stars have high /, it might explain why 
such modes have low amplitudes which are difficult to detect from 
the ground but are se en in space observations of these stars (e.g. 
iDegroote etal]|2009d) . 

In the past, model calculations for /? Cep and SPB stars have 
been mostly confined to modes of low spherical harmonic degree, 
/ ^ 2. With the precision attainable by space observations, it is 
expected that modes of much higher values of I should be de- 
tected. One may therefore expect that the view of /? Cep/SPB pul- 
sations from space may differ considerably from the view obtained 
by ground-based observations. This is indeed the case, as will be 
shown. 

The iron-opacity peak also leads to the formation of a small 
convective zone in the envelope of sufficiently luminous OB stars 
(iStothers & Chin 1993). It is often accompanied by an even smaller 
convective zone due to partial ionization of helium. These con- 
vective zones are so thin that they have no influence on the pul- 
sational stability or energy transport. The convection zone due to 
the iron-opacity bump is more prominent for lowe r surface grav- 
ity, hi gher luminosity and higher initial metallicity (Cantiello et al.l 
I2OO9I) . These thin convective zones may be connected to photo- 
spheric turbulence (as seen, for example in the broad spectral lines 
in O-type stars), spectral line variability, wind clumping and the 
generation of a magnetic field. 

In this context, one of the most interesting and important re- 
sults is that obtained from CoRoT observations of the /3 Cep star 
HD 180642 (V1449 Aql). In addition to a bulk of coherent modes 
with stable amplitude s, broad structures in t he frequency range 9- 
22 d^^ were detected jBelkacem et al .l2009h . These structures have 
been suggested to be modes with short lifetimes which are stochas- 
tically excited in a thin convective envelope. In addition, CoRoT 
observations of the O-type star HD46149 reveal a pattern charac- 
teristic of stochastic oscillations, which are predicted to be simi- 
larly ex cited as in B stars du e to thin convective layers in the en- 
velope ( IDegroote et al.l 120101) . With the outstanding precision of 
Kepler observations, it is possible to search for the presence of 
stochastically-excited modes in a larger sample of early B-type 
stars. 

Finally, of course, there is the prospect for uncovering hith- 
erto unsuspected types of variability among B stars. While it is not 
possible to deduce the nature of stars solely from an analysis of 
the light curves, one may obtain some clues by classifying stars 
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into similar groups from which a common mechanism may, with 
further observations, be deduced. In this paper we confi rm some 
of the results obtained by CoRoT dPegroote et al ] |2009ah . though 
testing of our speculations must await further observations. 

In this paper we review the photometric variability of all Ke- 
pler Asteroseismic Science Consortium (KASC) targets that were 
identified as main sequence B-type stars. The paper is organized as 
follows. Section [2] describes the selection of B-type stars from the 
Kepler photometry database. Section |3] shows how stellar parame- 
ters were derived for the selected sample. Section |4]describes the 
time-series analysis followed by a presentation and review of the 
different types of variability among the selected sample (Section[5]l. 
In Section[6]we discuss the SPB and /3 Cep variables that have been 
identified and compare their properties with models. The issue of 
stochastically-excited modes in early B-type stars is discussed in 
Section|7] 



2 SELECTION OF MAIN-SEQUENCE B-TYPE STARS 

Main sequence B-type stars (MSB stars) are members of Popula- 
tion I which are confined to the Galactic disk. Consequently, very 
few such stars are expected to be located at distances larger than 
about 100 pc above and below the Galactic plane. The Kepler field 
of view lies within Galactic latitudes 6° < 6 < 20° . A mid-B star 
with an apparent magnitude of 1^ = 10 mag at these latitudes will 
be about 600 pc above the Galactic plane at a distance of about 
4 kpc. Clearly, only the nearest and brightest B stars are expected 
to be included in the Kepler field of view. Only some of the bright- 
est stars in this field have been observed spectroscopically, which 
makes it very difficult to select the MSB stars. 

Fortunately, in preparation for the mission, Sloan Digital Sky 
Survey (SDSS)-like griz, D51 (510 ni n) and 2MASS JH K pho- 
tometry was obtained for most stars dBatalha et al.ll2010l) . Effec- 
tive temperatures, surface gravities and radii derived from this 

ahotometry are included in the Kepler Input Catalogue (KIC) 
. It should be noted that the above filters do not include the 
Balmer-jump band which is essential for the characterization 
of hot objects. Since effective temperatures are crucial for se- 
lecting B-type stars, we plotted the effective temperatures for 
stars with known spectral types to test the validity of the KIC 
calibration. The results shown in Fig. [T] indicate that the KIC 
temperatures are not very reliab le for B stars and are usually 
largely underestimate d (see also iMolenda-Zakowicz et alj l2010l : 
iLehmann et d] I2OI0I) . The recently released GALEX database 
( |http ://gal ex.stsci.edu/GR6/) can also be useful in 
distinguishing between B-type and later stars, although there is, as 
yet, few matches with B stars in the Kepler field (see also Sec- 

In order to maximize our chances of finding main sequence 
B-type stars in the Kepler database, we used the following selec- 
tion criteria: (i) all stars classified as B type in the HD or other 
catalogues and (ii) all stars with effective temperatures greater than 
10000 K in the KIC. This was supplemented by recent spectro- 
scopic investigations of st ars in the Kepler field, such as those of 
iLehmann et alj ( 1201 Oh and lCatanzaro et al] ( l201Clh . These selection 
criteria will inevitably lead to the inclusion of hot compact stars, 
mainly hot subdwarfs and white dwarfs. We therefore used the list 
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Figure 1. The effective temperatures from the KIC catalogue as a function 
of spectral type. The line is a recent effective temperature - spectral type 
calibration iMamaiekll2O10l) . 
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Figure 2. The UCAC3 proper motions of a selected sample of presumably 
MSB stars (filled and open circles). For comparison, the location of ex- 
cluded hot subdwarfs (squares) and white dwarfs (diamonds) is shown. The 
eleven MSB stars which are likely members of Gould Belt are shown with 
open circles. Four stars that lie outside the figure are shown with appro- 
priate symbols. The lai'ge circle with a radius of 6 mas yr^^ is centered at 
{{fia COS <5)o, figfi) = (-4.3,-5.4) masyr-i. 



of compact objects compiled bv lOstensen et alj ( I2OI0I) and elimi- 
nated all objects that were classified by these authors as non-main 
sequence stars. Several obje cts were excluded bec ause their effec- 
tive temperature derived by ILehmann et al.l ( I2OI0I) appeared to be 
smaller than 10000 K. These criteria resulted in a list of 72 pre- 
sumably MSB stars. 

Proper motions can be used to further refine our sample 
of stars. We used pr oper motions from the UCAC3 catalogue 
dZacharias et al.ll2009l) to identify and exclude probable non-MSB 
stars. All but seven stars in our sample are listed in this catalogue. 
Fig.|2]shows the location of the selected stars in the proper-motion 
diagram. As can be seen from the figure the MSB stars, espe- 
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Figure 3. The A/,i parameter as a function of {g — i) colour index for a sam- 
ple of MSB stars, white dwarfs and subdwarfs from the list of 0stensen et 
al. (2010). The symbols have the same meaning as in Fig.|2] Large symbols 
denote MSB stars with known MK spectral types. The dashed line, with the 
adopted ratio Ag/E{g — i) = 2.222 is drawn to separate MSB stars from 
compact stars. 

cially the brightest stars (with small error bars), form a well-defined 
clump centered at ((/.ic cos 5)o, /ii,o) ~ (—4.3,-5.4) masyr^^. 
The center of the clump is consistent with the average proper mo- 
tion of the young population in the direction of the Kepler field 
of view as can be judge d from proper motions of open clusters 
jKharchenko et alj|2003t) . However, there is a sample of eleven 
bright stars, shown with open circles in Fig.|2] which fall well out- 
side the clump but are classified as main sequence B8 or B9 stars. 
Their magnitudes and proper motions are consistent with member- 
ship of th e Gould Belt wh ich happens to cross the Kepler field of 
view (see lSlias et al.ll2009l and references therein). 

The three groups of hot stars we are considering here (MSB 
stars, hot subdwarfs and white dwarfs) show different spreads in 
proper motions (Fig.|2)- In particular, white dwarfs have the largest 
proper motions. The proper motions of hot subdwarfs are smaller, 
but still larger than for MSB stars. However, proper motions alone 
are not sufficient to distinguish between the tliree groups. For ex- 
ample, some hot subdwarfs are to be found among the clump of 
MSB stars in Fig.|2] Statistically, the larger the proper motion, the 
nearer the star. We may therefore use proper motion, ^, as an ap- 
proximate indicator of the distance, D. From this, we may define 
a quantity which is related to the absolute magnitude. This is use- 
ful because it is the absolute magnitude which best discriminates 
between the three groups of stars. 

We may define a parameter as an indicator of the absolute 
magnitude of a star: 

Mf. = 3 + 5 log(AAt/(mas yr-^)), 

where <? is a Sloan-like g magnitude taken from the KIC catalogue 
and A/i is the difference in proper motion from the center of the 
clump in Fig.|2] In order to avoid singularity in calculating Af^, the 
lower limit for was set to 1 masyr"^. The definition of is 
a slight modification of the parameter called "reduced proper mo- 



tion" bv l0stensen et alj ( 1201 Ol) for a similar purpose. A plot of 
as a function of the [g — i) index from the KIC is shown in Fig.[3] It 
can be seen that the three groups are quite well separated by A/^j . In 
particular, the dashed line in Fig.|3]separates MSB stars (above the 
line) from hot subdwarfs (below the line). The separation is verified 
by the location of MSB stars with good MK spectral classification 
which we know are certainly not compact objects (larger symbols 
in Fig.[3]l. Stars belonging to Gould Belt have different proper mo- 
tions and are located far from the clump of MSB stars in Fig. |2] It 
is therefore not surprising that they might have a larger A/^j and in- 
deed, KIC 1 1973705 which is a known MSB star (see TableQ, lies 
on the fringe of the Gould Belt distribution and below the dashed 
line. The spectrum of KIC 12258330, another star below the line, 
shows that it is an MSB star (Table[T). 

All stars falling below the dashed line in Fig.[3]except for these 
two stars were excluded from our sample as they are unlikely to be 
MSB stars. In addition, we excluded seven very faint stars with g 
magnitudes in the range 14.5-17.0 mag which do not have proper 
motions in the UCAC3 catalogue. It is unlikely they are MSB stars 
as they would lie at a very great distance from the Galactic plane. 
We are then left with a sample of 48 MSB stars for detailed fre- 
quency analysis. These stars are listed in Table[T] We have to stress, 
however, that even though the selection of our target stars was done 
in a consistent way, we cannot exclude the possibility that some 
MSB stars might have been overlooked. 



3 DETERMINATION OF FUNDAMENTAL 
PARAMETERS 

Stellar parameters are required in order to place the stars of Table[T] 
in the theoretical HR diagram. We cannot use the KIC parameters 
for reasons already discussed, but fortunately more reliable values 
of the effective temperature, Tch, and suiface gravity, logg, may 
be derived from spectroscopy (Section|3TT). Narrow-band photom- 
etry can also be used to derive the effective temperatures and lumi- 
nosities (Section [T2t . In addition, spectroscopic observations allow 
estimates of the projected rotational velocity, v sin i, which can be 
used as a constraint in the inteipretation of Kepler light curves. 



3.1 Spectroscopy 

In columns 2, 3 and 4 of Tabled we give Tefi, log g and v sin i, re- 
spectively, as derived from spectroscopic f ollow-up observations . 
The entries include the pub lished results of iLehmarm et al.l ( I2OI0I) 
and lCatanzaro et al] ( 120101) . New spectral classifications for some 
objects were also obtained with the B&C spectrograph using a 
400/mm grating (R ~ 550, A = 362-690 nm) attached to the Bok 
telescope at Kitt Peak observatory. These observations were done 
on June 14—16 and July 28-30, 2008. The Tefi and log g were de- 
rived by fitting model atmosphere grids to the hydrogen and helium 
lines visible in the spectra as described by 0stensen et al. (201(1). 
M etal-line blanketed LTE models of solar composition as described 
in iHeber et al.l ( I2OO0I) were used. Formal fitting eiTors on Tctt and 
log g are about 200 K and 0.05 dex respectively, while the system- 
atic en^ors are probably about 5 percent for both parameters . Be- 
cause of the low spectral resolution, the projected rotational veloc- 
ity could not be determined and v sin i = was assumed during 
the fitting procedure. 

For five targets (KIC 3865742, 6954726, 8057661, 8087269 
and 10536147), at least one high-resolution spectrum (i? = 81 000, 
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Table 1. List of 48 main sequence B-type stars discussed in this paper. The duration of the short-cadence (SC) and long-cadence (LC) runs for each star is 
given. The SDSS g — i index and the A/p parameter is discussed in the text. The notation FG means a star with frequency groupings while Rot implies the 
variability could be due to rotation. 
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Classincation 


Sp. type 


Notes 


number 


number 


[dj 


[dJ 


[mag] 


[mag] 


[mag] 








3240411 


19310672-1-3819417 


9.7 


322.1 


10.08 








Hybrid 


B2V'' 




3756031 


19350581-t-3852500 





321.6 


9.86 


-0.47 


10.44 


SPB + (5 Set? 


B5 V-IV * 




3839930 


19130682-1-3859204 




321.6 


10.68 


-0.51 


12.82 


SPB 


_ 




3848385 


19233702-1-3859363 




44.4 


8.86 


-0.35 


10.65 


Rot + ? 


B8V<= 


GALEX mags. 


3865742 


19411999+3856436 


_ 


321.6 


10.96 


-0.43 


11.03 


FG 


B 




4276892 


19400775-1-3920382 


25.9 


44.4 


9.06 


-0.35 


11.89 


Bin + ? 


B8-A0'* 


SB 1 ? FUV GALEX mag. 


4581434 


19455619-1-3939557 


30.0 


33.4 


9.12 


-0.20 


14.90 


Bin/Rot 


B9-A5 


SBl** 


5130305 


19563429-1-4014532 


_ 


321.6 


10.07 


-0.36 


13.51 


Rot 


B9 IV-V * 




5217845 


19555093-1-4023268 





321.6 


9.35 


-0.27 


12.32 


Bin + Rot 


68.5111* 


binary 


5304891 


19551942-1-4029573 


30.0 


44.4 


9.07 


-0.31 


12.98 


Rot 






5458880 


19390756-1-4037365 


28.9 


311.1 


7.63 


-0.53 


9.05 


Bin/Rot 


BO III "=,60.5 11 




5479821 


19572477-1-4038084 





321.6 


9.74 


-0.42 


12.63 


Bin + Rot/? 


B5.5V'' 




5786771 


19214208-1-4105482 


30.0 


44.4 


9.01 


-0.41 


14.13 


Rot 


B9-A5 




6128830 


19334968-1-4128452 


28.9 


44.4 


9.09 


-0.46 


12.19 


Rot 


B6HgMn'' 




6954726 


19344680-1-4226081 





321.6 


11.72 


-0.52 


13.62 


FG 




Be star, StHal66 


7599132 


19155121-^4312084 


33.4 


321.6 


9.22 


-0.44 


15.08 


Bin 


B8.5V'' 




7974841 


19463656-H4345480 


27.1 


44.4 


8.05 


-0.41 


14.88 


Rot/SPB 


B8V' 




8018827 


19115652-H4353258 


33.4 





7.90 


-0.45 


12.52 


Rot 


B9^ 


FUV GALEX mag. 


8057661 


20001799-H4349558 





321.6 


11.62 


-0.05 


12.80 


Hybrid 


OB-H'" 


Be star, KW97 41-38 


8087269 


191I0975-H4359571 


30.0 


321.6 


11.94 


-0.47 


14.99 


FG 




GALEX mags. 


8161798 


19213515-H4403022 





321.6 


10.33 


-0.44 


12.94 


Bin 


B8^ 


GALEX mags. 


8177087 


19422573-H4401183 


33.4 


311.1 


7.90 


-0.50 


11.16 


SPB 


B5V\B7ni 




8324268 


19565014-H4416159 





311.1 


7.81 


-0.42 


13.00 


Rot 


B9V^,B8pSi' 


V2095 Cyg, ct^ cVn star 


8351193 


19012177-1-4423115 


30.3 


33.4 


7.40 


-0.51 


12.73 


Rot 


B9V' 




8381949 


19484697-1-4421208 


_ 


321.6 


10.97 


-0.36 
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S ource of informatio n: KIC. ''iLehmann et alJj2O10l), ''ISato & Kuiill99 0). Cata nzaro et al.l bOlOl). ''iGuetteil lll968h. ^iHill & Lvnas-GravHl977|y 



g |Stephensonl il986h . ''lOownes & Keyed <1988l). 'iBartayj il983l). JlMacRaei jl952l) ''iHardorn et al. | 
" iBeckeJ h947l) . " iLindofj h972l) . p |0stensen et alj j2010l) . g |Svolopouloj il969t ) 



Il964. TGrenier et i 
. '' iDworetskyTTan " 



al.Nl999l) . ™ HP catalogue. 



Note: * Kepler magnitude. 



A = 370-1000 nm) was obtained with the Echelle SpectroPolari- 
metric Device for the Observation of Stars (ESPaDOnS) attached 
to the Canada-France-Hawaii Telescope (CFHT) at IVlauna Kea 
(Hawaii) on 2010 July 25, 3 1 and 2010 August 4. The Kepler light 
curves of the majority of these objects show characteristics similar 
to those of Be stars (cf. Section [5!3] l. Apart from the confirmed Be 
star KIC 6954726, none of these stars showed emission lines. These 



spectra were used only to derive a vsmi value by applyi ng the 
Fouri er transf orm method to the He I, IVIg I I and Si III lines ( Gray! 
ll982L see also lSimon-Diaz & Herrerol2007i for a recent application 
to OB-type stars). A full analysis of these spectra, together with 
a detailed analysis of their Kepler light curves, will be presented 
elsewhere (Gutierrez-Soto, in preparation). 
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Table 2. The determination of the stellar parameters (effective temperature T^ff, surface gravity logg, projected rotational velocity vs\ni and reddening 
E{B — V)) for the objects of Table[T]for which ground-based follow-up observations are available. Errors in these pai'ameters are given in brackets. We give 
the results obtained from spectroscopy in columns 3-5, from Stromgren photometry in columns 6-7 and from spectral energy distribution fitting in columns 
8-9. For each object, the log(L/Lg) given in column 9 was derived from the leftmost values of T|,ff and log g. 



Spectroscopy Stromgren photometry SED fitting 

KIC Tefi logg(cgs) usini log3(cgs) T^ff E{B-V) log(L/LQ) 

number [K] [dex] [kms"!] [K] [dex] [K] [mag] [dex] 
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Source of information: " iLehmann et alj j2010h . '' ICatanzaro et alj bOlOh . = Kitt Peak observations. CFHT observation. 



3.2 Photometry 

3.2.1 Stromgren photometry 

Standard Stromgren-Crawford uvbyP photometry was obtained 
with the 2.1-m telescope at McDonald Observatory. A stan- 
dard two-channel photoelectric photometer was used in 2010 Au- 
gust/September in single-channel mode. Apertures of 14.5 and 29" 
were used depending on the brightness of the target and sky back- 
ground. As the filter wheel can only hold four filters at once, uvby 



measurements were obtained at different times from /3 observa- 
tions. 

A set of standard stars, selected to span the whole parameter 
range of the targets in terms of (6 — y), mi, ci, /? and E{b — y), 
were observed for transforming the measurements into the stan- 
dard system. Transformation equations, used in the context of a 
larger observing program, were adopted from Handler (2010, in 
preparation). To maximize the number of objects for whic h stellar 
parameters can be estimated, we applied the calibration of ISalonj 
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j 19941) to the de-reddened Stromgren data to obtain Tcff and log g 
(see columns 5 and 6 in Table|2l(. As can be seen from Table|2] the 
results are, for the most part, in good agreement with the spectro- 
scopic values. 

3.2.2 Spectra energy distribution (SED) fitting 

Values for Teff were determined by fitting iKurucj ( 1 19931) model 
atmospheres with solar metallicity to the combined set of SDSS, 
2MASS and Stromgren colours. The calibration zero points for 
Stromg ren and SDSS were deri ved from the dedicated Vega model 
from iCastelli & Kurucj J 1994) and the 2MASS zero points were 
taken from Cohen et al.l 1 20031) . Whe rever possible , the p h otom- 
etry w as e xtended with va l ues fr om iDroege et al.l l l2006l) . lOfekl 
( |2008|) and lThompson et alj ( Il978l) . A grid search was performed, 
within the effective temperature range 8 000-30 000 K in steps 
of 200 K and the surface gravity range between 3.5 and 4.5 dex 
(cgs), in steps of 0.25 dex. For each combination of Tcff and logg, 
the model atmospher e was reddened according to the reddening 
law oflCardelli et al] ( Il989h using Rv =3.1. The colour excess, 
E{B — V), was determined by an iterative procedure using rnin- 
imization. The confidence interval for Tefi was chosen to encom- 
pass all values for which was below twice the minimal value. 
Equal weights were given to all data points in the minimization 
procedure. The value for log g cannot be obtained by this method. 

We tried to make additional use of the near- and far- 
ultr aviolet photometry from the Galaxy Evolution Explorer satel- 
lite jMartin et a"l]|2005h . The GALEX satellite was designed to im- 
age light from faint moderate- and high-galactic latitude objects 
in both near- and far-UV (NUV and FUV, respectively) imaging 
cameras by means of a beam-splitter. The latest General Release 
(GR6) delivery to the Multi-Mission Archives at Space Telescope 
Science Institute (MAST0 in 2010 showed an overlap of some 
60 percent of the Kepler field of view (FOV), which is far more 
than was delivered in previous releases. The GALEX GR6 cata- 
log were cross-matched with KIC stars at MAST for coordinate 
separations of less than 2.5 arcsec on the sky, resulting in FUV 
and/or NUV magnitudes for ten targets in Table [T] The SED fitting 
was redone including these GALEX magnitudes but did not lead 
to better results. Although GALEX-determined colours hold out 
great promise for differentiating among hot stars in the Kepler field, 
so far only explorations have been carried out to calibrate spectral 
types l lBia nchi et al..2007i) . Unfortunately, these are insufficient for 
our purposes. 



3.3 Luminosity 

From the definition of surface gravity, and eliminating the radius 
using the luminosity and effective temperature, we have 

log(L/Lg) = log(M/M0) - logg + 41ogTeff - 10.605. 

High-dispersion spectroscopy gives values of log Tcff and log g, 
but we need log(M/Mg) in this formula to calculate the luminos- 
ity. To do this, we used the relationshi p expressing l og(M /Mg ) in 
terms of log Tcff, logp and [Fe/H] in [Torres et al.l 1 I2OIOI) (setting 
[Fe/H] = 0.0). Since the values of Tcff and log g as derived from 



^ Support for the MAST is provided by NASA through a grant from the 
Space Telescope Science Institute, which is operated by the Association of 
Universities for Research in Astronomy, Inc., under NASA Contract NAS5- 
26555. 
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Figure 4. The theoretical j3 Cep and SPB instability strips (from 
iMiglio et ai]l2007h showing the location of Kepler targets for which the 
stellar parameters were derived from ground-based follow-up observations 
(Table|2j. Large filled circles are Kepler SPB//9 Cep hybrids. Open cir- 
cles are Kepler SPB stai's and the asterisks ai'e stars with frequency group- 
ings. The small filled circles are non-pulsating Kepler stars and the cross is 
KIC 1 1973705. The squares with filled circles are two CP stars. The filled 
square is VI 449 Aql. The location of the empirical & Scuti instability strip 
as determined from Kepler S Set stars is indicated by the contour in the bot- 
tom right comer The error bars at the bottom left-hand comer indicate a 
typical standard deviation of 300 K in effective temperature and 0.5 mag in 
absolute magnitude. 



the spectra are probably more accurate than those obtained from 
the Stromgren/Balona calibration and SED fitting, the log(L/Lg) 
values were derived from the spectroscopic determinations where 
possibl e. Failing this, values determined from uvbyP photometry 
and the lBalonal7l994l) calibration were used. The adopted values 
of log(L/Lg) are given in the last column of Table |2] As men- 
tioned before, the spectroscopic and photometric values are in good 
agreement. 

Comparison of the resulting positions in the HR-diagram with 
known /3 Cep and SPB stars and with the instability strips for modes 
of low and high degree is shown in Fig.|4] 



4 DATA AND FREQUENCY ANALYSIS 

Following each quarter of a 370-d Kepler solar orbit, the tele- 
scope is rolled to keep its solar panels facing the Sun. Choices and 
changes of targets are made on the basis of these "quarters". Most 
data are obtained in long cadence (LC) with integration times of 
29.4 min. For 512 targets data are obtained in short cadence (SC) 
mode with sampling ti mes just under 1-min . Characteristics of LC 
mode are d escribed inljenkins et al.l ( 1201 Obh and characteristics of 
SC mode in lGillilandetal.l ( l2010h . An overview of the K epler sci- 
ence processing pipeline is given bv ljenkins et al. I( l2010ah . The du- 
ration of the quarters are as follows: 
QO: 10 days (2009 May 2-2009 May 11), 
Ql: 33 days (2009 May 13-2009 June 15), 
Q2: 89 days (2009 June 19-2009 September 16), 
Q3: 89 days (2009 September 18-2009 December 16), 
Q4: 90 days (2009 December 19-2009 March 19). 

These observations provide a unique opportunity to study the 
incidence of pulsation at the lowest amplitude levels. The median 
noise level, (Tm, in the periodogram depends on the stellar magni- 
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10 

Frequency [d"'] 

Figure 5. Illustration of the determination of the false alarm probability 
(FAP) for KIC 1 1 360704. The periodogram based on the Q2-Q4 data before 
and after prewhitening of 40 frequencies is shown in the top and bottom 
panel, respectively. 



tude, V , and on the duration of the mn, At. We find that a relation 
of the form: 

log(crM/ppm) = 0.2(5 - 8.0) - 0.5 log(At/d) + 0.9, 

where ctm describes the median noise level of the Kepler peri- 
odograms quite well. A more practical measure of detection thresh- 
old, however, is the level of the top-of-the-peaks in the periodogram 
(the grass level), which is approximately equal to (3-4)ctm. For our 
sample of B stars and full Q0-Q4 coverage this detection threshold 
is typically equal to about 1 ppm for the brightest stars and about 
15 ppm for the faintest. 

In order to determine whether a peak in the periodogram is sig- 
nificant requires the calculation of a false alarm probability (FAP) 
which is related to the periodogram noise level. A determination of 
the FAP, st r ictly v alid only for equally-spaced data, was obtained 
by IScargld dlgsl . Scargle's FAP requires the existence of a set 
of frequencies at which the periodogram powers are statistically 
independent. The existence of these is guaranteed only when the 
data are evenly spaced in time. Although the Kepler data are al- 
most consistently evenly spaced, the time series is interrupted for 
scheduled spacecraft rolls. Further, observations of individual stars 
will contain odd outlier points that need to be eliminated before 
calculating the periodogram. Inspection of quarters Q2-Q4 of the 
Kepler data (see below) show that approximately 5 percent of the 
evenly-spaced data set is lost in this way. The conditions required 
for Scargle's FAP are therefore not strictly met. 

Even if the Scargle conditions had been met, Scargle's method 
would require us to calculate the periodogram power only at the 
natural frequencies. This restricts candidate pulsation frequencies 
to the Scargle grid. The resolution of this grid is much poorer than 
that allowed by the length of the Kepler data set, sacrificing preci- 
sion in the pulsation spectrum. Preliminary determinations of the 
pulsation spectra of the stars discussed in this paper show a large 



number of very closely spaced frequencies. The restriction imposed 
by t he Scargle method is th us excessively debilitating in this case. 

iFrescura et al. ] ( l2008h showed that when the FAP is deter- 
mined by Monte Carlo methods, independent frequencies are not 
needed, allowing the periodogram to be sampled on arbitrary fre- 
quency grids. They also showed that oversampling of the peri- 
odogram in the frequency domain has no deleterious effect on fre- 
quency determination. We compu ted FAPs for Ke pler quarters Q2- 
Q4, using the methods detailed in lFrescura et al.l ( i2008.) . The FAPs 
determined by this method are shown in the bottom panel of Fig.[5] 
where they have been superimposed on a heavily prewhitened pe- 
riodogram (with 40 frequencies) of KIC 11360704 for illustrative 
purposes. The non-prewhitened periodogram for this star is shown 
in the top panel of Fig.|5] In this figure, the two lowest FAP levels 
lie so close to the frequency axis in this plot that they can not be 
usefully displayed. The FAP levels can therefore be more clearly 
seen when plotted in Fig[5] lower panel. This shows the extremely 
high statistical significance of practically all the visible peaks in 
the top panel of Fig. [5] and of many dozens of peaks still present 
in the lower panel of Fig. |5] A detailed analysis of the statistical 
properties of the Kepler data is in preparation and will be published 
separately. 

The data used here are from the short commissioning run (QO) 
and the first four quarters (Q1-Q4). Data for some stars are avail- 
able for all these runs, a total observing duration of over 320 d in 
long cadence. Periodograms using the long-cadence data have far 
better resolution and are more useful in determining the frequencies 
and amplitudes. SC observations were mostly used to search for 
frequencies above 24.5 d~ ^ , the Nyquist frequency for the LC data. 
We did not find frequencies higher than the long-cadence Nyquist 
frequency in any star except for KIC 11973705. The duration of 
the observing run for each star in short- or long-cadence mode is 
indicated in columns 3 and 4 of Table[T] 

The MAST Kepler archival database provides for corrected 
and uncorrected data. We used both corrected and uncorrected data, 
depending on a star and the observing quarter. Before the data were 
subjected to an analysis of the frequencies, the time sequence was 
examined and corrected for trends and discontinuities in the light 
curve. To a large extent this is an automated procedure in which 
a smooth curve is fitted to the time series and residuals examined 
for systematic departures indicative of discontinuities. At the same 
time, points which are clearly outliers are rejected. While this auto- 
mated procedure worked satisfactorily in most cases, the data cor- 
rected in this way were visually examined and the corrections done 
by hand if deemed unsatisfactory. Owing to various instrumental 
effects, such as aperture changes, sudden pointing deviations, long- 
term drifts and many other effects, the data contain signals which 
do not belong to the star. Some of the signals lead to characteristic 
peaks in the periodogram which have been tabulated by the Kepler 
team and are well-known. Nevertheless, the behaviour at low fre- 
quencies is still poorly understood. In analyzing these data, we take 
the view that frequencies below about 0.2 d~^ should be regarded 
with suspicion. 

It should be noted that while these corrections are often sub- 
jective in nature, they only affect the low frequency range. Fre- 
quencies above 0.5 d~^ are hardly affected at all by the corrections 
and comparison of frequencies done independently by different au- 
thors are in good agreement. One might also question the level at 
which a certain peak is deemed to be sig nificant. In som e analyses 
the Lomb-Scargle false alarm criterion (IScarglelll982l) was set to 
some level, above which peaks are regarded not to be significant. 
In others, the local signal-to-noise level was calculated and the peak 
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Figure 6. Fourier frequency spectra of Kepler LC data for seven SPB//3 Cephei hybrid stars. Two periodograms are presented for each star: calculated for 
original data (left) and after prewhitening with a number (indicated) of low-frequency terms (middle). The periodograms were calculated up to the Nyquist 
frequency (24.5 d^^) for the LC data. The panels on the light show histograms of weighted frequency spacings (HFS) plotted as a function of spacing, A/. 
The numbers in a spacing bin are nomialized to the highest bin. 



deemed significant only if it exceeds a certain factor. These deci- 
sions do not affect the periodogram, but only the set of frequencies 
which are thought to be real. In this paper, the frequencies we will 
be discussing are certainly significant by any criterion. 



Characteristic frequency spacings in the periodogram, if they 
exist, clearly provide information about the star. In order to detect 
characteristic frequency spacings we calculated histograms of all 
possible frequency separations between unique frequency pairs af- 
ter the frequencies were extracted. It is necessary to impose some 
frequency and amplitude restrictions in using this method. This was 
done by attributing arbitrary weights to the spacings. The weights 
depend on the amplitude, Wij = \og{Ai ■ Aj), where Ai and Aj 
denote amplitudes of the two terms that form a spacing and Wij 
is the weight attributed to the spacing A/ij = |/i — A bin 
of 0.01 d~^ was used in these histograms which will be hereafter 
abbreviated as HFS. 



5 RESULTS OF THE ANALYSIS 

The periodograms and light curves of the stars in Table [T| were ex- 
amined and deductions made which are described below. It is, of 
course, difficult to deduce the true nature of the star just from the 
periodogram. Nevertheless, with this information supplemented by 
the stellar parameters that were derived for most stars under inves- 
tigation, we have made an attempt at understanding the true nature 
of detected variations. 



5. 1 SPB//3 Cephei hybrids 

When we examined the periodograms of Kepler B stars we were 
immediately struck by their strange appearance. In ground-based 
observations, and even observations from MOST and CoRoT, a 
P Cep star is easily recognizable in having distinctive peaks in the 
high-frequency domain. It has long been assumed that the typical 
frequencies for j3 Cep stars are in the range 4-20 d~^ with perhaps 
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Figure 7. Top: The same as in Fig.[6]but for SPB stars. Bottom: Frequency spectra for three stars with a small number of detected terms. 



a few low-amplitude peaks in the region 0.5-2 d^^ which charac- 
terize typical P Cep/SPB hybrids. The SPB stars, on the other hand 
show distinct peaks in the low-frequency domain and no frequen- 
cies much higher than 2 or 3 d~^. Here we assume that stars with 
frequencies exceeding 3.5 d~^ can be regarded as belonging to the 
(3 Cep group since this frequency is the approximate boundary be- 
tween p and g modes in non-rotating stars on the zero-age main 
sequence. We realize this is an arbitrary number, but we chose it 
simply for classification purposes. 

The periodograms of Kepler pulsating B stars are completely 
different. They show a very large number of peaks in the low- 
frequency domain, but with some isolated low-amplitude high- 
frequency peaks. Nothing like this has ever been seen before and 
we need to ask ourselves how such stars should be classified. They 
are not typical /3 Cep stars because the scattered high-frequency 
peaks usually have very low amplitudes and the periodogram is 
dominated by a forest of low-frequency peaks. They are not typi- 
cal SPB stars either because SPB stars do not have low-amplitude 
high-frequency peaks. They do not resemble any of the known hy- 
brids which are really (3 Cep stars with a few isolated, weak, low- 
frequency peaks. A possible solution to the strange appearance of 
the periodograms is discussed in Section |6] 

There are seven stars in our sample that we will call 
SPB//3 Cephei hybrids to stress that the dominant modes are the 
SPB-like g modes (Fig.[6ll compared to the jS Cephei/SPB hybrids 
where the jS Cephei p modes have the highest amplitudes. The peri- 
odograms of the two types of hybrids are very different. Since low 
frequencies dominate in all seven Kepler stars, the middle panels 



of Fig. |6] show frequency spectra prewhitened with the strongest 
low-frequency modes in order to reveal the high frequency peaks 
more clearly. One assumes that the low frequencies are g modes 
while the few isolated high frequencies may be normal p modes. 
However, this assumption may not be correct because the effect of 
rotation may scatter some low-frequency g modes into the high fre- 
quency domain (as measured by the observer). For example a sec- 
torial mode of low intrinsic frequency having / = 20 and azimuthal 
number m = 20 may well appear at a frequency of 10-20 d^ ^ since 
the rotational frequency perturbation is proportional to mfi where 
S7, the rotational frequency, could be around 1 d^^ in some stars. 

The detected amplitudes are very small and except for the 
dominant period in KIC 10960750 do not exceed 700 ppm. As was 
indicated above, however, the 320-day long Kepler data are of un- 
precedented precision in which the noise level is only several ppm. 
The number of significant modes that can be extracted is there- 
fore very large and greatly exceeds the number of frequencies in 
ground-based observations of hybrid stars. The histograms of fre- 
quency spacings are shown in the right-hand panels of Fig.|6] These 
histograms basically reflect the structure of the periodogram in a 
different way and are quite interesting for tracing similarities be- 
tween one star and the next. It is interesting to note the presence 
of a 2/ harmonic of the prominent high-frequency peak at / = 
10.449 d"^ in KIC 10960750. 

Of the seven SPB/^ Cep hybrids, KIC 324041 1, 8057661, 
10536147 and 10960750, have known values of vsini (Table|2]l. 
The first two stars have low projected rotational velocities typical 
of SPB stars, but the last two stars have high v sin i. There are very 
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few SPB stars known from ground-based observations in which 
Dsini > 100 km s~^. In spite of the pronounced difference in 
projected rotational velocities between these two pairs of stars, we 
can find no obvious difference in either the appearance of the peri- 
odograms or their associated fr equency spacing distributions. 

KIC 8057661 (KW41-38, [Kohoutek & Wehmevej[T99^) . an- 
other SPB//3 Cep hybrid, is included in the catalogue of objects 
with Ha emission jH ardorp etaLlll964 . However, recent spectra 
do not show any emission in the Balmer lines and the spectral lines 
are relatively narrow (J. Gutierrez-Soto, in preparation). It may be 
therefore a Be star seen at low inclination which temporarily does 
not show Ha emission. 



5.2 SPB stars 

The next group of eight stars we discuss here are stars that can be 
called 'normal' SPB stars because they show only low-frequencies, 
being devoid of the isolated high frequencies in the stars discussed 
above. Their frequency spectra are shown in Fig. |7] The three fre- 
quency spectra shown at the bottom of Fig.|7]are for stars where the 
number of detected modes is too small to calculate a reliable HPS. 

Only three stars from this group have known v sin i values 
which are low in all three cases: 31 kms"^ for KIC 3756031, 
22 kms"^ for 8177087 and 53 kms"^ for 8459899. These values 
are quite typical for ground-based SPB stars. 

In two stars, KIC 8177087 and 10797526, the formal solution 
comprises of a large number of low-frequency terms. These two 
stars are, however, among the brightest in our sample. For all other 
bright stars (g < 8.5 mag) a similar increase of amplitude towards 
low frequency is observed. We therefore conclude that this increase 
of the noise towards low frequencies is due to the instrumental ef- 
fects caused by the fact that the stars are close to the bright limit of 
Kepler observations. 

There are several high-frequency terms in the frequency spec- 
trum of KIC 3756031. They are barely visible in Fig. |7] because 
their signal-to-noise does not exceed 4.5 (amplitudes 8-16 ppm), 
though this is still significant. The star is located in the middle of 
the SPB instability strip where p mode instability is not expected. 

There is another star which may perhaps be regarded as an 
SPB candidate, KIC 1 1973705. The frequency spectrum of the LC 
data (upper panel of Fig. [8} clearly shows multimode behaviour in 
two frequency regions: below 10 d~^ and over 18 d~^ with only 
some small-amplitude terms in between. There are SC data avail- 
able from QO which, however, lead to much poorer frequency res- 
olution. In pailicular, the group of modes around 3 d^^ is not well 
resolved in the SC data. Nevertheless, the data are useful to show 
that high frequencies extend beyond the Nyquist frequency of the 
LC data (/n = 24.5 d"^). For example, the peak at /{ = 20.68 d"^ 
in the LC data (Fig. [8} is not real; the true frequency of this mode 
is /i = 28.26 d^^ as we see clearly from the SC data. The LC 
frequency is mirrored with respect to /n : /i — /n ~ /n — .f'l ■ 

Th e star was found to b e a double-lined spectroscopic bi- 
nary bv iLehmann et al.l ( |2010|) . The variability of KIC 11973705 
is dominated by an almost sinusoidal variation with period P — 
3.3857 d and semi-amplitude of about 3560 ppm (Fig.llOt. This is 
very probably the orbital period or half this period. The high fre- 
quencies probably belong to the A-type star (Tcff — 8 000-9 000 K, 
^ehmann et al. 2010) which is then a S Set variable. One could in- 
terpret the low frequencies as SPB pulsations of the B-type primary, 
but since many 5 Set stars have low frequencies (7 Dor-type vari- 
ability), these could equally well belong to the A-type companion. 

The frequency spectrum of KIC 10658302 shows weak evi- 



dence for two broad peaks, at about 1.5 and 3 d~^. It may perhaps 
be classified as one of the stars with characteristic frequency group- 
ings discussed below. 



5.3 Stars with frequency groupings 

There is one ce rtain Be st ar in o u r sample, KIC 6 9 54726 
= StHa 166 jStephensoiJ 1 19861 ; iDownes & Kevej 1 19881 ; 
I Gutierrez- Soto et al.l I2OI0I) . Its projected rotational velocity 
(160 km s~^. Table |2l( is moderate, so that it is probably a rather 
low-inclination Be star. Its periodogram (top panels of Fig. [9]l 
shows a high-amplitude peak at 1.027 d^^ and three broad peaks 
(0.1, 0.9 and 1.7 d^^). The HFS for this star shows maxima for 
separations A/ ~ 0.2, 0.8 and 1.5 d^^. In fact, the light curve 
of this star is dominated by quasiperiodic long-term variability 
which was subtracted at the preprocessing phase. The frequency 
spectrum for this star seems to be characteristic of Be stars. 
For instance, very similar groupin gs were discovered in several 
Be st a rs observed by the M OST Iwalker et af] |2005| ; ISaio et al 



2007"; 'Cameron et al." 200 8|) and Coieor satellites jHuat et al 



2009j ; ,Diago et al. 2009 ; Neiner et al .1120091; ICutierrez-Soto et al 
20ld) . These groupings are usually interpreted in terms of high 



radial order g modes that have frequencie s in the corotating 
frame smaller than the rotational frequency ( Walker et al.l |2005| ; 
iDziembowski et al] 1 20071 ; ICameron et al. Another expla- 

nation is that these short-term variations are due to surface or 
circumstellar inhomogeneities (lBalonall2009h where the second 
and third broad peaks are interpreted as the rotational frequency 
and its first harmonic. A combination of both effects is also a 
possibility. In both hypotheses, the separation of the broad peaks is 
interpreted as the rotational frequency, fi. 

Seven other stars from our sample show frequency spectra 
which, to some extent, are similar to that of KIC 6954726, i.e., 
they show at least three frequency groupings as judged from their 
HFS (Fig.|9]l. None of these stars is known as a Be star. In the HR 
diagram (Fig.|4l, they cover the range of effective temperature cor- 
responding to B2-B8. Projected rotational velocities are available 
for four of these stars and are moderate to high, ranging between 
133 and 271 km s"^ (Table|2l(. For comparison, the mean projected 
velocities of B6-B9 main sequence stars ranges from 73 km s~ ^ 
(B9-9.5 IV) to 144 km s'^ (B6-B8V) jLevato & Grossoll2004l) . 
The critical rotational velocity for these mid- to late-B stars is about 
400 km s"\ The mean equatorial velocities of Be stars is about 
350 km s~^. Our sample of these stars is too small to draw any re- 
liable conclusion, but it would be indeed interesting to verify if the 
rotational velocities of stars in this group are systematically higher 
than the mean for B-type stars. 

If either of the hypotheses proposed for the short-term vari- 
ability in Be stars is also valid for stars with period groupings, the 
spacing between frequency groupings. A/max, can be interpreted 
as fi. Thus it may be possible to estimate the rotational frequency 
independently of spectroscopic observations of v sin i. As a result, 
the inclination of the rotation axis, i, can be derived. In a pulsating 
star, knowledge of i might help in mode identification and calcula- 
tion of the visibility of the modes. The values of A/max from the 
HFS for all eight stars shown in Fig. |9] are presented in Table |3] 
We note that for two stars, KIC 8766405 and 12258330, two differ- 
ent values of A /max were derived. The former shows a very com- 
plicated HFS consisting of broad peaks with superimposed much 
narrower ones, the HFS for the latter has a double peak. 

We suspect that the origin of the frequency grouping in these 
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Figure 8. Frequency spectra for Kepler LC (top) and SC (bottom) data of KIC 1 1973705. The highest peak at 0.2953 d is truncated. The vertical dashed 
line corresponds to frequency 24.5 d~^, the Nyquist frequency for the LC data. 



Table 3. A/max derived from the HFS for eight Kepler B-type targets. 



KIC 




KIC 




number 


[d-i] 


number 


[d-i] 


3865742 


0.78 


9716301 


0.24 


6954726 


0.79 


10285114 


0.52 and 2.76 


8087269 


1.59 


11360704 


2.05 


8766405 


0.88 


12258330 


2.24 and 2.95 



stars could be rotation, though the sample at our disposal is too 
small to verify this idea. 

5.4 Stars with variations due to binarity or rotation 

About a half of stars in our sample (23) show much less compli- 
cated light curves. Their variability can be understood in terms of 
proximity effects in a binary system, the presence of surface in- 
homogeneities and rotation or the combination of both. With the 
photometric precision of the Kepler data, one may expect to de- 
tect weak proximity effects such as l ight reflection, def ormation 
of components and Doppler beaming jZucker et al.ll2007l) . even in 
low-inclination systems. The light curves of these stars are shown 
in Figs.[TO]and[TT] The periods and amplitudes of the main period- 
icity (A\) and its second harmonic {A2) are given in Table|4] Here 
we have included KIC 1 1973705, the binary in which the compan- 
ion is a 5 Scuti variable. 

None of the light curves of stars in this group shows evi- 
dence of eclipses. A thorough discussion and analysis of the bi- 
nary light curves is beyond the scope of this paper and would 
probably require time-series radial velocity observations. Nonethe- 
less, four stars from this group (KIC 4276892, 4581434, 5217845 
and 12207099) are known (or suspected) to be binaries (Table [T]!. 
Without corresponding radial velocity observations we cannot be 
sure that the photometric period corresponds to the orbital pe- 
riod. We also suspect that other stars showing very regular vari- 
ability with relatively large amplitudes, especially those that have 
double-wave light variations, are probably binaries as well. These 
include KIC 5458880, 7599132 and 8161798 (Fig. [Hi and also 
KIC 5479821 and 9716456 (Fig.QIll. 

The longest period detected in our sample is 12.9 d in 
KIC 4276892, a suspected SBl binary l lCatanzaro et"ani2010l) . Bi- 
naries with longer orbital periods may have a photometric signature 
in the Kepler data, but if the amplitude is too small it will probably 
be mistaken for an instrumental trend. 



There are several variables (KIC 4581434, 5130305, 5786771, 
8018827, 8351 193) with periods in the range of 0. 4-0.6 d. Among 
these, only KIC 4581434 is known to be a binary jCatanzaro et al.l 
I2OIOI) . All but KIC 8018827 show variability with the peak-to-peak 
amplitude of less than 30 ppm. This very small amplitude may be 
due star spots or other inhomogeneities in a rotating star. If the 
star spots migrate with respect to each other this could explain am- 
plitude modulation in the light curves. Amplitude modulation is 
present in the five stars mentioned above. Their frequency spec- 
tra (Fig. not show residual power in the vicinity of the subtracted 
main periodicity and its harmonics. KIC 8018827 is exceptional be- 
cause the strongest residual signal is observed at a frequency of 
2.269 d^^, i.e., about 10 percent less than the dominant frequency. 

Similar variability can also be seen in stars that have longer 
periods, e.g., KIC 5304891 and 8389948. For stars where i; sin i is 
known, an additional argument in favour of rotation as the cause 
of variability is the fact that the observed periods are consistent 
with vsini. With the assumption that the main periodicity coixe- 
sponds to the rotation period, the minimum radius can be estimated: 
R-mi-a = Protvsm ill-K. The values of i?min are given in Table |4] 
for all stars that have measured vsmi. We note that the radii of 
late B-type stars on the main sequence range from about 1.5 to 
3.5 R0. On this basis, rotation as the cause of variability seems to 
be excluded for KIC 1 1973705, the longer period of KIC 5479821 
and also perhaps for KIC 5217845. The first star and the last star 
are known binaries, so we suspect that their variability originates 
from this cause. 

KIC 5479821 is one of two stars (Fig. II 111 which shows two 
non-sinusoidal periodicities. If the longest period in these two stars 
is attributed to orbital effects, the shortest period needs to be inter- 
preted differently. Rotation is one possibility, but if this is the case, 
the systems are strongly non-synchronously rotating. Contamina- 
tion or multiplicity cannot be excluded either. Variability caused 
by rotation and binarity are not mutually exclusive. In short-period 
binaries, synchronization of the orbital and rotation period can be 
expected, but because B-type stars are relatively young, some devi- 
ations from full synchronization is pos sible. In this context, w e note 
that KIC 5217845, which is a binary jLehmann etal]|2010l) . has a 
light curve which is strongly modulated with a period of about 6 
days. 

Among the stars discussed in this section there are two chem- 
ically peculiar (CP) stars. KIC 8 324268 (HP 189160, V2095 Cyg), 
was discovered as a CP Si star bv lZirinl(ll95lh . Its photometric vari- 
ability was detected in Hipparcos data, but the period could not be 
derived. In Kepler data the star shows strictly periodic variability 
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Figure 9. The same as in Fig.|6]but for eight stai's showing frequency groupings. 



with a period of almost exactly 2 days (Fig. [TO] Table |4j and only 
very small amplitude modulation of about 0.2 percent. The period 
is assumed to be the rotational period as in other CVn stars. 

The other CP star is KIC 6128830, a HgMn star 
dCatanzaro et al.l l2010l) . The light curve is nearly sinusoidal 
with a period of about 4.84 d and peak-to-peak amplitude of 
3600 ppm. The photometric period, if assumed to be the rotational 
period, is consistent with vsini —15 kms~^ (Table [4l>, whic h 
is typical for the slowly rotating HgMn stars dAbt et al]|l972h . 
Light variations in HgMn stars have on l y rece ntly been detected 
in CoRoT observations dAlecian et al ] l2009h and line profile 
variations are know n only in very few stars (iHubrig et al.l 1200^ ; 
iBriguet et alj|201ol) . The light variation implies an inhomogeneous 



distribution of various elements in the atmospheres of HgMn 
stars which challenges our understanding of the nature of these 
stars. HgMn stars were known not to have strong magnetic 
fields l lAuriere et af] 1201 Ol) and, while diffusion should lead to 
vertical stratification, no horizontal surface abundance patches 
are expected. However, recent observations dHubrigetal.lbOld) 
indicate the presence of a weak longitudinal magnetic field in the 
HgMn star AR Aur 

As shown in Fig. |4l stars from this group are mostly late 
B-type stars falling outside both the P Cephei and SPB instabil- 
ity strips. This may explain the lack of variability other than that 
related to rotation or binarity in these stars. Those that are lo- 
cated within the instability strips are either CP stars or do not 
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Figure 10. Phase diagrams (left panels) and frequency spectra of the residuals (right panels) after removing the dominant periodicity shown in the left panel. 
In the phase diagrams the means in 0.01 phase intervals are also plotted. The long ticks at the top of frequency spectra indicate the frequency of the principal 
periodicity and all its detected harmonics below 10 d~^. The ordinate in all panels is given in ppm. 



show variations that can be interpreted as due to pulsations. Only 
KIC 5479821, which is located in the SPB instability strip, shows a 
single peak at 5.709 d^^ and amplitude of 8 ppm which is difficult 
to understand. The variability in KIC 11973705 has already been 
discussed. 



On the other hand, there are frequencies that cannot be eas- 
ily explained in four late B-type stars: KIC 3848385 (0.31 d"^), 
4276892 (1.29 and 1.45 d"^), 7974841 (0.28 and 0.52 d"^), and 
8692626 (6.09 d'^). 

Finally, we are left with two stars, KIC 9655433 and 



Kepler observations of variability in B-type stars 15 



30000 
10000 


-10000 

5000 





-5000 



p -iiai32a6i d 

A A A 

; \ i \ i \ 



ilOOO 
20OO 


-3000 
HiO 

Am 

200 


-200 

-Am 




P - 5^130956 d 



P - 1.11523^1(1 



■ ■■■■■ 



0.0 0.5 OjO 0.5 OjO 0.5 0.0 0.5 

Phase mase 



Ml 'I I' 
1 11 11 



KIC5i|79a21 
OO- CM(LjC) 



KIC971&I56 
OO - CM(LjC) 



23i|567a9 10 
Pret^enoy [d"'] 



Figure 11. The same as in Fig.[T()] but for two stars, KIC 5479821 and 9716456, with two non-sinusoidal periodicities detected in their light curves. 



Table 4. List of the parameters of 24 B-type stars with vaiiations caused primarily by binarity and rotation. 





Period 






V sin % 


^' mm 


Notes 


number 


[d] 


[ppm] 


[ppm] 


[kms-i] 


[R©] 




3848385 


0.9689(7) 


9.5(7) 


3.1(7) 








4276892 


12.85859(6) 


581.4(5) 


1194.6(5) 


10 


2.53 


SBl? 


4581434 


0.502901(8) 


19.2(5) 


11.6(5) 


200 


1.98 


SBl 


5130305 


0.464508(7) 


17.5(4) 


4.8(4) 


155 


1.41 




5217845 


0.839114(2) 


5643(9) 


38(9) 


237 


3.91 


binary 


5304891 


0.98750(3) 


34.5(14) 


37.7 (14) 


180 


3.49 




5458880 


3.51193(4) 


2710(16) 


4675(16) 








5479821 


1.7011778(4) 


13998.4(19) 


3429.5(19) 


85 


2.84 






5.430956(19) 


3043.2(19) 


86.1(19) 




9.07 




5786771 


0.403879(6) 


28.1(5) 


3.6(5) 


200 


1.59 




6128830 


4.84129(3) 


1688.5(18) 


205.9(18) 


15 


1.43 


CP 


7599132 


1.303685(3) 


8890(20) 


8458(20) 


63/50 


1.61/1.28 




7974841 


2.8426(5) 


107.9(26) 


41.1(26) 








8018827 


0.3978805(9) 


437.3(4) 


415.2(4) 








8161798 


2.2029818(4) 


13208(4) 


39551(4) 








8324268 


2.0091295(4) 


13185.7(7) 


990.0(7) 






CP 


8351193 


0.56767(4) 


1.90(24) 


2.10(24) 








8389948 


1.00602(3) 


12.98(26) 


4.50(26) 


142/130 


2.81/2.57 




8488717 


3.25485(7) 


287.7(12) 


45.2(12) 








8692626 


1.6483(3) 


4.86(27) 


1.23(27) 








9716456 


1.8132861(7) 


859.0(11) 


4455.4(11) 










1.115234(25) 


84.6(11) 


17.3(11) 








10130954 


2.246541(10) 


406.7(5) 


106.5(5) 








11973705 


3.38570(8) 


3558(14) 


33(14) 


103 


6.85 


SB2 


12207099 


1.42256744(14) 


16159.4(7) 


3283(7) 


43 


1.20 


SB2? 


12217324 


6.557(7) 


53.3(8) 


4.1(8) 









1 1817929, which we have not yet discussed. The former is the only 
one from our sample which shows no evidence of variability in Ke- 
pler data. The latter is a member of the open cluster NGC 681 1 and 
shows the frequency spectrum characteristic of 5 Set stars. Its spec- 
tral type is very poorly known (Table [T) and we suspect that it may 
be an A-type star. 



6 MODES OF HIGH DEGREE IN B-TYPE STARS 

We have seen that there are several B stars in the Kepler data base 
which are clearly SPB or SPB//3 Cep hybrids according to our defi- 
nition. We have already remarked on the vastly different appearance 
of the periodograms which are so different from those in ground- 
based observations of these stars. This requires an explanation. 

We see from Fig. |4] that all of the SPB//3 Cep hybrid stars are 
within the ji Cep instability strip, with KIC 8714886 just cooler 



than the red edge, but within the SPB instability strip. The instabil- 
ity strips generally discussed in the context of ground-based obser- 
vations are constructed for modes with i ^ 3. From the ground we 
only see modes of high amplitude, and these tend to be modes with 
/ ^ 2 because cancellation effects render modes with higher I invis- 
ible in photometry. The amplitude limit for ground-based observa- 
tions is typically a few millimagnitudes, but for Kepler photometry 
the detection limit is only a few ppm. At this level, the visible range 
of / must be considerably increased. Cancellation effects are severe 
up to Z ~ 4, but f or higher values of I the decreas e in visibility with 
/ is much smaller jBalona & Dziembowskil 19991) . It is possible that 
with Kepler data, pulsations with values of I as high as 20, or even 
higher, may be visible. Why a particular mode should be selected in 
preference to another is, at present, an unsolved problem and there 
is no reason to suppose that only modes of low degree are present 
in these stars. Clearly, we cannot limit the comparison with models 
to low values of /. 
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One of the puzzling aspects of the Kepler SPB//3 Cep 
hybrids is the fact that their frequency spectra are dominated 
by low frequencies. On the other hand, ground-based obser- 
vations of SPB//3 Cep hybrids are dominated by the high 
frequencies with only a few low-frequency SPB-type modes. 
iBalona & Dziembowskil ( 1 19991) have studied the excitation of 
modes of high I in B-type stars. They find that for a given value of I 
there are two separate frequency ranges where unstable modes may 
occur: a range of SPB-like low frequencies and a high-frequency 
(/3 Cep like) range. In the high-frequency range only modes of low 
I are excited, while in the low-frequency range modes of high I are 
excited. These high-degree modes are trapped in the outer layers. 

In zero-age main sequence stars the region of SPB-type low- 
frequency instability is present for stars with 4. 1 < log(Tcft /K) < 
4.4 (B7 to B I) and reaches its maximum extent for Z = 8. As the star 
evolves, the SPB-type region extents further to hotter temperatures 
up to log Tefi = 4.5 (09, BO), reaching its maximum extent for / = 
10. The low-frequency region is still present, though smaller, for 
I = 24. It turns out that the maximum extent of the instability strip 
for low-frequency g modes for / up to 24 more or less coincides 
with the SPB instability strip for low / shown in Fig.|4] It therefore 
covers almost the entire B star region. Thus it is not surprising to 
find that most of the stars observed by Kepler show many low- 
frequency modes. In all likelihood, these are g modes with high /. 
They are not seen from the ground because of their low amplitudes. 

It would seem that none of the stars observed by Kepler are 
in the region of the HR diagram where p modes are predominantly 
excited. The hottest stars in our sample all lie on the cool side of the 
p-mode /3 Cep instability strip as can be seen in Fig.|4] The high- 
frequency modes are perhaps only weakly driven and have very 
low amplitudes. Of course, rotation also plays a big role, especially 
for modes of high I which can be perturbed to higher and lower 
frequencies in proportion to the rotational frequency. The peculiar 
appearance of the periodograms of these stars thus has a logical 
explanation in terms of excitation of high-degree modes and their 
location on the cool side of the /3 Cep instability strip. 

In Fig.|4] we also show the location of Kepler B stars which 
do not pulsate. As already mentioned, some of these stars may 
be binaries, but there is nothing that clearly distinguishes pulsat- 
ing and non-pulsating stars. In ground-based observations, the SPB 
stars nearly all have wsini < 100 km s^^, which is not a selec- 
tion effect and suggests that high rotation suppresses pulsation. The 
Kepler sample of pulsating and non-pulsating stars comprise both 
slow and rapid rotators, so rotation does not appear to be a factor. 
We note that some of these non-pulsating B stars are within the 
SPB strip and two are within or close to the j3 Cep instability strip. 
This indicates that the instability strips are not pure, i.e. that stars 
need not pulsate even if they are located within the instability strip. 



7 STOCHASTIC EXCITATION OF MODES 

As we mentioned in the introduction, one of the most inter- 
esting results from space observations of /3 Cep stars was the 
possible discovery of completely unexpected modes with short 
lifetimes i n CoRoT observations of the B1.5II-III /3 Cep star 
yi449Aql (Belkace m et al.l2009l) and in the O-type star HD 46149 
jPegroote e t al. 201(J. Features in the periodograms of VI449 Aql 
were found in the range 9-22 d^^ and were interpreted as stochas- 
tic modes excited by turbulent convection j ust below the photo- 
sphere in the thin Fe a nd He convective zones dCantiello et alj2009l ; 
iBelkacem et al ] [2O10l) . The occurrence of convection in the region 



of iron opacity bu mp depends strongly o n the luminosity and metal- 
licity of the star dCantiello et"ai]|2009l) . One thus expects that the 
occurrence of stochastic oscillations in B stars would follow a sim- 
ilar dependence if the modes are excited in the Fe-bump region. 

Confirmation of the existence of stochastically excited modes 
in massive stars is clearly of great interest. One of the problems 
facing such an investigation is that it is not a simple matter to dis- 
tinguish between self-excited and stochastic modes, as well as non- 
pulsational light variations. If oscillations occurred in the same fre- 
quency range as the /3 Cep pulsations, as they do in V1449 Aql, then 
one cannot immediately eliminate the Fe bump «: mechanism as the 
excitation mechanism for these modes. The suspected stochastic 
modes seen in the Coi?or observations appear as wide structures in 
the periodograms of V1449 Aql and HD 46149. They can either be 
interpreted as a large number of self-excited modes driven by the 
n mechanism (possibly beating with a rotation period unresolved 
in the measurements) or as single modes that change amplitude and 
phase stochastically. It is not possible to distinguish between these 
two ex planations unless further information is available. 

In lBeikacem et alj feoogi) the evidence is presented as a time- 
frequency diagram in which the constant amplitude of a nearby 
self-excited mode is contrasted with the variable amplitude of the 
modes thought to be stochastically excited. These authors find that 
the stochastic modes are organized as regularly spaced patterns in 
the periodogram, which is stated to be signature of these modes. 
Regular spacing applies only in the asymptotic limit of high fre- 
quencies and occurs for both stochastic and self-excited modes. The 
claimed modes in V 1449 Aql may well be stochastically excited, 
but unfortunately it cannot be proven at this stage. What lends cre- 
dence to such an hypothesis, however, is that there is a prediction 
of where stochastic modes are exp ected to occur in the HR diagram 
among the B stars jCantiello et al.,.2009.) . It would certainly be im- 
portant to search for wide structures in the periodograms of the 
Kepler B stars to determine if such stars conform to this prediction. 

We visually examined the periodogram of CoRoT observa- 
tions of V 1449 Aql to familiarize ourselves with the broad struc- 
tures. We then looked for similar structures in the periodograms of 
the pulsating stars shown in Fig. |6l but nothing resembling them 
could be found. All the stars have well-resolved peaks in the re- 
gion 10-24 d^^. Moreover, the background noise level conforms 
to the relationship discussed in Section|4]except for KIC 3865742 
and KIC 10960750 where it is about a factor of 3 or 4 higher than 
expected for the stellar magnitude and length of the time series. 
For both stars this is a result of the huge number of low frequen- 
cies. There are at least 600 signals with frequencies less than 6 d~^ 
in the Kepler light curves of these two stars. 

KIC 5458880 has an effective temperature and luminosity 
very similar to V 14 49 Aql (Tefi = 24 500 K, log 3 = 3.45, 
iMorel & AertsI I2OO7I) . It has a strange double- wave light curve 
which we have interpreted as a proximity effect of a binary. The 
periodogram shows a noise level which increases towards low fre- 
quencies with no clear evidence of distinct frequencies, except at 
the binaiy period and its harmonic, down to a level of about 10 ppm. 
There is certainly no sign of broad structures. The peaks we exam- 
ined have amplitudes of about 10 ppm, i.e., an order of magnitude 
smaller than those found in VI 449 Aql. 

Our conclusion is that the existence of stochastic modes in B 
stars is not yet proven and that the Kepler stars show no indication 
of the presence of these modes. If the modes do exist, this result 
may simply be due to the fact that the Kepler stars in which /? Cep 
pulsations are seen are too cool to have a sufficiently active con- 
vection zone. The occurrence of stochastic oscillations relies on the 
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occurrence of convection in the iron opacity bump region. Since we 
cannot yet explain tlie occurrence of constant stars within the jS Cep 
instability strip, it is also possible that a star may not pulsate even 
if it is in a location in the HR diagram where stochastic driving is 
expected to be efficient. 



8 CONCLUSION 

Kepler observations of 48 main sequence B stars show that 15 
of them pulsate with frequencies typical of SPB stars. Of these 
stars, seven show isolated, weak, high frequencies in the /3 Cep 
range. The frequency spectra differ considerably from those seen 
in ground-based SPB//? Cep hybrids and can be understood if most 
of the low-frequency modes are of high spherical harmonic degree. 
Models show that low-frequency modes of high spherical harmonic 
degree, I, up to Z = 24 are excited in stars spanning practically the 
whole main sequence band in B stars. Since we expect to see far 
more modes of high / in Kepler observations than in ground-based 
observations, this may explain why so many low-frequency modes 
are visible in the Kepler data. Although the presence of low- and 
high-frequency modes place these stars in the same class as the 
ground-based SPB//3 Cep hybrids, it is not clear whether this clas- 
sification is meaningful in the context of space observations. We 
call the hybrids observed by Kepler SPB//3 Cep hybrids to empha- 
size the dominance of modes in the SPB region. 

We also find that though nearly all 48 B stars in our sample 
are variable at some level, there are certainly stars within the SPB 
as well as the cooler side of the /3 Cep instability strip which are 
not SPB stars by the usual definition of these stars. This indicates 
that there are mechanisms which restrict pulsations in B stars. Such 
mechanisms may be rotation (which appears to inhibit pulsations in 
SPB stars), magnetic fields and chemical stratification. It would be 
important to determine which (if any) of these factors are responsi- 
ble for damping of pulsations by a comparative study of pulsating 
and non-pulsating stars. We may also conclude from these observa- 
tions that the j3 Cep instability strip is not pure, i.e., there are stars 
in the fi Cep instability strip that do not pulsate at all (KIC 5458880 
and 10130954). This conclusion is strengthened by the B0.5 IV star 
HD 51756 for which the variations observed in the CoRoT light 
curve are interpreted as due to rotational modulation (Papics et al., 
submitted). There are also many stars within the SPB instability 
strip which do not seem to pulsate. 

We have not found any indication of pulsating star situated 
between the hot end of the 5 Set instability strip and the cool end 
of the SPB instability strip. This conclusion is not only a result of 
the study of B stars reported here, but also follows from a study 
of main sequence A s tars in the Kepler field not reported here. 
lOegroote et al.l ( l2009al) found a number of such stars. One star that 
show clear high-amplitude 5 Set pulsations, KIC 11973705, is al- 
most certainly a B star with a 5 Set companion. 

We have identified a group of stars which show a structure 
in the periodogram consisting of unresolved frequency groupings 
which may or may not be due to pulsation. The light curve shows 
a period corresponding to the mean frequency of the first broad 
peak (not counting the lowest frequencies) with a strongly modu- 
lated amplitude. The structure of the periodogram in many cases 
resembles that of Be stars, but those stars with spectra show that 
they are not Be stars. Moreover, a similar structure is found among 
Kepler A and F stars as well. The HPS is a valuable tool in studying 
such structures and may possibly be applied in other cases as well. 
However, the physical meaning of these structures is not presently 



understood, and without further study it does not seem worthwhile 
to pursue this line of investigation. 

We have looked for the broad structures seen in CoRoT ob- 
servations of the /3 Cep star V1449 Aql and the young 0-type 
binary H D 46149. They were in terpr eted as solar-like stoch astic 
modes by iBelkacem et al.l | |2009|) and iBelkacem et alj | |2010| ), re- 
spectively. We are unable to find such structures in the Kepler stars, 
but this may be because they are too cool. A problem is that it 
is not really possible to distinguish between quasi-periodic (non- 
pulsational) variations, stochastic oscillations and closely-spaced 
unresolved self-excited oscillations. A conclusive demonstration 
that stochastic modes are present in B stars is not possible in prin- 
ciple with existing data. However, it is possible to test the hypothe- 
sis that convectively-driven stochastic modes occur within a certain 
region in the HR diagram by comparing the location of stars with 
broad features in the HR diagram to the predictions. If the compar- 
ison supports the hypothesis, then the concept of stochastic pulsa- 
tions in B stars would certainly become more acceptable. At this 
stage, such a comparison is not possible with Kepler data. 

It will be very difficult to use the observed frequencies on their 
own for asteroseismological modelling without mode identification 
for at least some frequencies. Because of the low amplitudes, mode 
identification requires multicolour space observations which will 
not be available in the foreseeable future. Some progress might be 
possible if systems of equally spaced frequencies or periods can 
be found. A preliminary investigation suggests that this might be 
possible in one or two stars. If this can be confirmed, it might tell 
us something about the rotational splitting or perhaps even lead to 
an estimate of the large separation. 

Finally, we have detected variability in many stars which we 
suspect may be proximity effects due to binarity. In other stars, the 
modulated light curves are more easily explained in terms of mi- 
grating star spots. We have also detected periodic light variations in 
two chemically peculiar stars. Kepler observations of B-type stars 
have produced some suiprises as well as new opportunities for a 
better understanding of these stars. 
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